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Abstract 

A measurement consists in coupling a system to a probe and reading the output of the probe to 
gather information about the system. The weaker the coupling, the smaller the back-action on the 
system, but also the less information conveyed. If the system undergoes a second measurement, 
the statistics of the first output can be conditioned on the value of the second one. This procedure 
is known as postselection. A postselected weak measurement of an observable can give a large 
average output of the probe when the postselected state is nearly orthogonal to the initial state of 
the system. This large value is an interference effect in the readout of the probe, which is initially 
in a coherent superposition of readout states (also known as pointer states). Usually, the weak 
interaction between system and probe is considered instantaneous, so that the dynamics of the 
probe can be neglected. However, for a weak measurement in solid-state devices, an interaction of 
finite duration is likely needed. 

Here we show how this finite duration generates a contribution of the dynamical phase to the 
readout statistics. Furthermore, we derive interpolation formulas that are able to describe the 
statistics of the weak measurement for the whole range of pre- and postselected states. Phase- 
space averages appear in the expansion, suggesting an interpretation in terms of non-positive 
probabilities. Decoherence in the probe is also accounted for and it is pointed out the existence of 
a regime of intermediate coupling strength in which coherent oscillations can be observed in the 
probability of the readout. 
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I. INTRODUCTION. 



Sequential or joint nonprojective measurements can reveal the quantum nature of the 
detectors. For example, the pioneering work of Arthurs and Kelly [1] demonstrated how 
two weak measurements of momentum and position could be carried out simultaneously, 
and how the uncertainty principle obeyed by the detectors (or probes) showed up in their 
mutual back-action. While the theory of positive operator-valued (POV) measures (of which 
the weak and the strong measurements are particular cases) has its mathematical roots in 
the seminal papers of Neumark [2, 3] and has been developed since the 70s [4-6], only in 1988 
was it realized that a weak measurement followed by a strong one could lead to arbitrarily 
large values of the average output [7], provided that the weak measurement is conditioned 
on the result of the strong one. This conditioning is called postselection. As postselected 
weak measurements can give an arbitrarily large average output (even for non-orthogonal 
preparation and postselection [8]), they have been used to amplify a weak signal [9-11], to 
settle fundamental issues — such as determining the traversal time of a barrier [12], solving 
Hardy's paradox [13, 14], and observing quantum trajectories [15]- - and also to perform 
quantum state tomography [16] . Concerning the last application, there are several proposals 
| l] for extending the procedure to mixed states. It has also been shown [21, 1 !] that 
having a stream of particles sent to probes initially prepared to measure in the strong regime 
can create coherence in the probes; this drives the measurement to the weak regime and 
is reflected in a deviation of the statistics from that predicted for projective measurement, 
provided the decoherence rate of the probes does not exceed the firing rate of the particles. 
For other applications and studies of the weak measurement, see Ref. [23]. 

The theoretical papers on postselected weak measurements are mostly limited to the study 
of the average value and assume an instantaneous (von Neumann) interaction. While for 
optical implementations postulating an instantaneous interaction is reasonable, for the still 
prospective realizations in solid-state systems (see the proposals [24-30]) it is more realistic 
to suppose a coupling that lasts a finite time. In an earlier paper [31], we introduced a finite- 
duration interaction, assuming only that the measurement was a quantum nondemolition 
one [32]. The results of Ref. [31], however, were limited to the average value of the output and 
its variance (barring the case study of a spin 1/2, for which, since an exact solution exists, 
the expansion in the weak coupling was made merely to test its validity). In a more recent 
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paper [33], one of us studied the whole statistics of the weak measurement, by performing a 
controlled expansion in the coupling strength, as in Ref. [34], and providing an interpolation 
formula that works for any preparation and postselection. Here, we provide the statistics 
of a nondemolition weak measurement of an arbitrary variable. We show that the finite 
duration has observable consequences, since in a weak measurement the coherence of the 
probe manifests through the contribution of its density matrix off-diagonal elements (in the 
readout basis), and since the interaction contributes dynamical phases to these elements. 



II. DESCRIPTION OF THE MEASUREMENT. 

Let us consider a quantum system prepared, at time tj, in a state psiii) = Pi (preselection) 
and a second quantum system, the probe, prepared, at time to, m a state pp(t ) = p . Let 
the Hamiltonian of the system be Hs = hfls, and that of the probe Hp. The system and the 
probe interact, at time T« > tj, to, through Hi nt = —h\g(t)xA, where A is an operator on the 
system's Hilbert space, x on the probe's, and g(t) is a function that vanishes outside a finite 
interval [T^Tj-], with J g(t)dt = 1. The interaction generates a time-evolution operator U 
that acts on the Hilbert space of the system and the probe, entangling the two. Let the 
operator k be the conjugate observable [? ] of x, [x, k) = i. It follows from Heisenberg's 
equations of motion, or even from Hamilton's equations for the corresponding classical case, 
that k is the observable of the probe that carries information about the measured quantity A. 



A, He 



0, otherwise 



As the interaction lasts a finite time, A has to be be conserved, i.e. 
the question would arise of what is being measured [? ]. We notice that, in order for the 
measurement to be nondemolition, k too must be conserved during the free evolution of the 
probe and change only due to the interaction with the observed system. Hence Hp = ftwp(k). 
In other words, we are considering not a von Neumann (instantaneous) weak measurement, 
but the more general nondemolition (finite-duration) weak measurement, having the former 
as a special limiting case. At time t > Tf a projective measurement of the observable k 
is made on the probe. As k is conserved by Hp, the value obtained will not depend on 
the time, as long as t > Tf. In a strong nondemolition measurement, it is furthermore 
assumed that the probe is prepared in a state such that the elements (k,a\UpoW \k,a) are 
negligible unless k ~ k(a), where the function k(a) maps the values of a to the readout. This 
hypothesis guarantees the fidelity of the measurement, as well as its repeatability. However, 



System 




no, keep k 




Probe 




Interaction 



yes, discard k 



FIG. 1. A schematic view of the measurement with pre- and postselection, the horizontal direction 
representing increasing time. A random number generator yields < p < 1 after each trial. If 
P < W (S) the outcome is considered. 

as we are interested in the weak regime, we drop this assumption. 

III. POSTSELECTION IN A MIXED STATE. 

At time tf > Tf, a projective measurement of an observable S of the system is made, 
giving an output S and leaving the system in the state \S). Then, given S, one keeps the 
outcome k according to an arbitrarily chosen probability w(S). This step leaves the system 
in the postselected mixed state 



The procedure detailed above and sketched in Fig. 1, describes a measurement with pre- and 
postselection. Without loss of generality, we shall consider ti = t = Tj = and tf — Tf — r, 
with r the duration of the interaction. We notice that W = Y2s w (S) is not necessarily 1 
as w(S) are probabilities conditional on the event S, not probabilities of the event. If the 
system has a finite-dimension Hilbert space, with D the dimension, then < W < D. For 
w(S') = wo if S = So and otherwise, the postselection is in the pure state p/ = |So) (Sol- In 
this case, w < 1 is a sub-optimal choice, in the sense that some of the trials are discarded 



Pf 



(1) 
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unnecessarily. The opposite limit is found for w(S) = w = W/D ,\/S. This implies that no 
postselection is made, i.e., pf = D~ L 1. The optimal choice is w — 1. If all the probabilities 
are multiplied by the same factor r, rw(S) —> w(S), the postselected state is unchanged, 
but the probability of a successful postselection is also multiplied by r, and so is the joint 
probability of observing the probe in p and successfully post-selecting the system. Thus the 
conditional probability, given by the ratio of the two probabilities above, is unaffected by this 
rescaling. The optimal choice of r is such that it maximizes the probability of postselection, 
namely max s [w(-S')] = 1. Hence W opt > 1, with the equality only for postselection in a pure 
state. In a sense, we could say that the probabilistic postselection leaves the system in the 
unnormalized state ^2w(S) \S) (S\. 

Another proposed method [35] is to make the postselecting measurement a POV one. 
This way p/ is replaced by a positive operator Ef. However, while Ef appears in the 
probabilities, as e.g. in Eq. (5) below, the system is not selected in a state Ef/TrEf but in 

A 1 /2 A 1 /2 

Pf = E/ PiEJ , with p its reduced density matrix after the interaction with the probe. 



IV. EXACT RESULTS. 



The joint state for the probe and the system, at any time t > r is 

p(a, a , k, k'; t) — J dkidk 2 po{ki, k 2 )pi(a, a) (k, a\ Wo,r \ki, a) (k 2 , a'\ 1A\ T \k', a) 

x exp(i{[w P (/0 - u) P (k)} + [w s (o') - u> s (a)]}(jt - r) ), (2) 



where U is the time evolution operator generated by H$ + Hp + H int , \k, a) are the simul- 
taneous eigenstates of k and A, and 0Js( a ) are the eigenvalues of Cl s corresponding to \a). 
Following a lemma demonstrated in the Supplemental Material section, there is an analytic 
solution for the propagator, 

(k, a\ W ,r \ko, ao) = S a<ao 5 (k - k - Xa) e - iVa{k \ (3) 

where we define the Hamiltonian phase 

r o (Jfe) := fdsu;p{k- \a[l-h(s)\), (4) 
Jo 

with h(s) = Jq ds'g(s'). The joint probability of observing the outcome k at time t for the 
probe and of postselecting the system in pf at time Tf, follows readily from Eqs. (2) and 



(3): 

r(k, p f ) =W ( a '\ Pf \ a ) ( a \ Pi K> e~ i[ra(fc) - r -' (fc)1 po(fc — Xa,k — Xa'). (5) 

a,a' 

For a von Neumann measurement, T a (k) = 0. 

The conditional probability of observing the outcome k, given that the state has been 
postselected in pf, follows from Bayes's rule, 

ew - (6) 

' post 

where the denominator V pos t = J dkV(k, pf) represents the probability of making a success- 
ful postselection in pf, irrespective of the value of k, or of what observable of the probe, if 
any, was measured. What remains to be done is simple: apply the controlled expansion of 
both V(k,pf) and V pos t in A as in Ref. [33], with the difference that here one should keep 
track of additional contributions from the dynamical phases T. It is fundamental to keep in 
mind that, both in V pos t and in V(k,pf), the zeroth and first order terms vanish for nearly 
orthogonal pre- and postselected states (NOPPS). Thus one should write down Q(k) as the 
ratio of two quadratic polynomials in A, without succumbing to the temptation to expand 
the denominator. 

V. DEFINITIONS. 

Before proceeding to the expansion, we introduce the normal weak values as in Ref. [33] 

a m>n = Tis{A m p f A n Pl }. (7) 

We note that a„ i?n = a* mn) so that a m ^ m are real. For NOPPS, ao,o — > and ao, n — > 0, 
with ao,o/tto,n ~ ^ for n > 1; also, a mtTl —¥ p ^ for m,n > 1. While there can be 
exceptions to this behavior when A n+l = A for some integer n, in any case ct^i stays finite 
for NOPPS, barring some trivial instances. This consideration is important, as a\ t i provides 
the dominant term to both V(k,pf) and V pos t for NOPPS. 

We shall also use the phase-space averages [36], f(x,k) = J dkdx f(x, k)H^ (x, k), with 
the initial Wigner function of the probe. In particular, we define the covariance C (/, g) = 
f9~f~9- -X natural concept arising in the expansion is that of phase-space conditional 
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averages f(x,k)> k , defined by 

P)|^oW = j dxf(x,k)U^(x,k) (8) 

with Vo(k) = po(k,k). The phase-space average of f(x,k) is found by integrating Eq. (8) 
over k. Finally, we define the time scales r n = f Q ds [h(s)] n \l — h(s)], which satisfy r n > 
T n+1 ,n E N. 

VI. WEAK COMPARED TO WHAT? 

As A is a dimensioned constant, it has to be compared to some homogeneous quantity 
in order to establish whether the measurement is weak, strong, or intermediate. Let K k be 
the coherence scale, i.e., k ^> n k ==>• Po(K + k, K — k) <C Po(K, K), and A k the classical 
uncertainty scale, i.e., k ^> A k ==>■ p (k,k) <C p (0,0). Precisely, we define K k = (A^.) -1 , 
where A x is the classical uncertainty of the conjugate variable of k, = x 2 — x 2 , and 
^■1 = k 2 — k . The uncertainty relation requires n k < 2A k . In order to realize the weak 
(coherent) regime, the coupling constant must be small compared to the coherence scale, 
precisely, Xcim Kk, with clm the maximum distance between the eigenvalues of A. The 
validity of the expansion relies also [34] on Ai < 1. Equation (5) shows that both the 
off-diagonal terms of the detector and the Hamiltonian phase contribute to the statistics. 
Thus we are in presence of interference. We also remark that if either pf oc 1 or pi oc 1, i.e., 
no postselection or no preselection is made, the coherent contributions disappear. Hence, 
both pre- and postselection are essential in order for interference to show up. 

VII. CONTROLLED EXPANSION. 

The probability of postselection is found by expanding the integral of Eq. (5) (for sim- 
plicity, W — 1), 

Ppost - a ,o - 2Ax ro Im(a ,i) + A^^ai,!. (9) 

where Xt = x + u' P (k)t and the prime stands for differentiation. Notice that u' P (k) is the 
velocity, so that Xt is the displacement of the x variable. Due to its non-trivial dynamics 
as k, and hence the velocity, changes during the interaction with the system, r appears 
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instead of r [? ]. Since the second-order term contributes significantly only for NOPPS, we 
have neglected a 0j 2 compared to a.\ tl . 

Furthermore, as all expressions are homogeneous in a m>n , we may reduce the number of 
independent parameters by dividing by, say, 0:0,0- Accordingly, we introduce the canonical 
complex weak value A w = a^i/c^o anc ^ the additional real term B w = o^i/ao.o- We 
remark that B w > \A W \ 2 . The equality holds whenever pf and pi are mixtures of pure states 
J2f w f \ f) (f\ an d J2i w i K) (*1 eacn P a i r °f which has the same weak value (/| A \i) / (f\ \i) = 
A w Vi, / : WiWf 7^ 0. In particular, B w = \A W \ 2 for pure pre- and postselected states. 

It is convenient to introduce the characteristic function Z{6) = J dke t0k Q(k), from which 
the moments of the distribution can be generated as (k n ) = {—id/d9) n Z(0)\$=q, with (• • • ) 
average over Q. By integrating Eq. (5) times exp[i6k], changing the variable of integration 
to k — X(a + o')/2, expanding, and normalizing, 



with 



NZ(6) ~ e idk + X i9e m Re(A w ) - 2e i6k x T0 lm(A w ) + \ 2 e iek (x* Q + jj B w , (10) 



N= I J!^i^ l-2\x T0 \m{A w ) + \ 2 xlB w . (11) 

"0,0 



The interpolating formula for the probability is instead 



NQ = V 



1 - 2\x T0 \ k Im(A w ) + \ 2 x 2 T0 ^ k B 



A 



V' Q Ke{A w )-^B w 



(12) 



The derivatives come from having expanded po(k — Xa, k — Xa'). This way, the shifts of the 
probability distribution are lost. However, as in the weak regime the shifted peaks are not 
resolved, this yields a small deviation from the exact result. 

As an example, away from NOPPS, we can expand the average value of k to first order, 

(k)~k + ARe(AJ - 2XC(x T0 ,k)\m(A w ). (13) 

The contribution from lm(A w ) is usually neglected, as in an instantaneous measurement 
x T0 = x and as a rule x and k are assumed to be initially uncorrelated. The correlator 
C(x, k) is proportional to the derivative of the variance of k under specific hypotheses [37]. 

Furthermore, it may be useful to observe the value x of the probe, instead of k. We give 
the full statistics in the Supplemental Materials section, and give here the average value of 
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k 

[- Kxact App'oximate \ 

FIG. 2. Comparison between the interpolating (dashed) and the exact formula (solid) for the 
probability Q{k) of obtaining k in a weak measurement of A = a z . The preselection is in a pure 
state oriented at an angle ir/3 with the z axis, and the postselection makes an angle it — 0.1 with 
the preselection and lies in the plane defined by m and a. All variables are in units of A^: A = 0.5, 
K/i = 2, ku = 10, kp = oo. A not so weak coupling strength was chosen, in order to have a 
discernible difference between the two curves. 

x, which is, in the regime away from NOPPS, where we are allowed to keep up to first order 
terms, 

(x) ~x T -2\ [C(x T ,x T0 )lm(A w ) - u/£t v Re(A w )] . (14) 
The results of Ref. [37] are recovered as a special case. 

VIII. DECOHERENCE. 

As an effective model for decoherence, we consider a random classical force acting on the 
probe, while the latter interacts with the system. In principle, if the probe is prepared before 
the interaction starts, decoherence will tend to kill the off-diagonal elements po{k, k'), driving 
the measurement to the weak incoherent regime. However, this effect can be counteracted 
by preparing the probe shortly before its interaction with the system, and in any case it 
can be treated rather simply. On the other hand, decoherence during the interaction with 
the system is unavoidable, as the probe is open to external influences, which can come from 
the system and from the environment. For simplicity, we take up(k) = hk 2 /2M. The effect 



-4 -2 2 4 

k 

[- Kxact Approximate \ 

FIG. 3. Coherent oscillations of the readout probability. Here, kjj = 0.2 and the other parameters 
are as in the previous figure. 

consists in the addition of an imaginary phase to T a — T a / in Eq. (5), which becomes 
V(k,p f )=^2(a'\p f \a) (a|p i |aOe^-(* ) - r -'^e- A2 ^ ') a7 * flT ^/ M po(A;-Aa,A;-Aa / ). (15) 

a,a' 

We assumed the random force to have the correlator (f(t)f(t')) = S(t — t^fikgT, with p a 
constant, T the temperature, and kg Boltzmann's constant, and we defined the decoherence 
rate 7 = p/M [? ]. The factor e is a dimensionless parameter that depends on the details 
of g(t). Thus there is a new scale for A to be compared with, the decoherence scale Kp = 
{fkBTr 3 /2M]~ 1 / 2 . If A ^> Kd, the contribution from the off-diagonal elements of po(k,k') 
becomes negligible, and we fall into the weak incoherent regime. If A Kr>, the net result 
of decoherence is simply to shift the coefficients of B w in Eqs. (10) and (12). 

IX. A CASE STUDY. 

We consider a spin-1/2 system preselected in a state pi = (1 + m ■ c)/2, on which 
a weak measurement of A = a • cr is made, and that is postselected in pf = (1 + n ■ 
<t)/2. The interaction is considered constant g(t) = r _1 and the probe is prepared in 
Po(k, k') oc exp { — (k + fc') 2 /8Aj* — (k — A; / ) 2 /2/t 2 }. The free Hamiltonian of the probe is as 
in the previous section h 2 k 2 /2M. This defines the Hamiltonian scale ks = \f2M/hr. We 
show the approximate and exact probability in Fig. 2. 
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In the intermediate case kjj/Ak <C A <C k^, it is no longer legitimate to expand exp (iT a ). 
Then, as shown in Fig. 3, the probability displays coherent oscillations in k. 

X. CONCLUSIONS. 

We have studied the statistics of a weak nondemolition measurement, providing expres- 
sions for the probability and the characteristic functions that are robust for any overlap 
between the postselection and the preparation, contrary to the results of Ref. [7] and sub- 
sequent papers on weak measurement. We have included decoherence in an effective, albeit 
heuristic, way, and we have pointed out the existence of a regime of intermediate strength, 
in which coherent oscillations can be observed. 
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Appendix A: A useful lemma 



We prove the following lemma: for a system subject to the time-dependent Hamiltonian 
H = hcu(k) — f(s)x, with u,f arbitrary functions and [x,k] = i, with the time-evolution 
being hence U = Texp < — i J^ds ui(k) — f(s)x \, the propagator reads 



(k\U\ko) = 5 ^k — k — J dsf(s)j exp —i J dscjUz — J ds'f(s' 



(Al) 



1. Brute force derivation 



We use the path-integral technique, but without path-integrals, i.e. we approximate the 
time-ordered exponential as a product of AT + 1 terms, and we introduce between each term 
the identity in the \k) basis, obtaining (&at+i = k) 



(k\U\ko) — / dki ■ ■ ■ dk,N J^J(^j+i| g x P "I — i£j w(k) — f(s)x 

f n 

~ / dk\ ■ ■ ■ dk^s exp [— ieju(k)} exp [iejf(tj)x]\kj) 

J j=0 

f n 
= / dki ■ ■ ■ dk N J^J 5(kj + i — kj — Ejf(tj)) exp [— i£jU)(kj + i)] 



3=0 
N 



8 j k - k - ^2 £ jf( t j) ) ex P \ ~i 
j=o 



N 



.3=0 



m=0 



(A2) 



In the limit A^ — > oo, the sums in the last line of Eq. (A2) become integrals, and the 
approximated equality with the first line becomes exact, thus the lemma is proved. 



2. Alternative derivation 

We provide an alternative, more elegant derivation of the above result. The technique 
illustrated below may find applications in other fields. The Schrodinger equation in wave 
number space is 

u(k)i>(k,t)+if(t)—i/>(k,t) = i-gt/>(k,t). (A3) 
13 



Upon rearranging the terms and dividing by ip(k,t), we have the non-homogeneous first 
order partial differential equation 



f(t) 



d__d_ 
dk~di 



u(k, t) = iu(k), 



(A4) 



with u = ln^>. The solution to the corresponding homogeneous equation is u yk + f Q f(s)ds 
with uq an arbitrary function of one variable. [? ] This suggests to change the variables to 
H,t, with n = k + J* f(s)ds. Let t>(/i, t) = u(k(/j,t),t). The PDE becomes then 



- J^G"' *) = i u ~ J f( s ) ds 



and a particular solution is readily found 

t) = —% J Loyt — J f(s')dsjds, 
so that, in terms of the original function, the general solution is 



(A5) 



(A6) 



u(k,t) = -i J ulk + J f{s')ds'jds + u \k + J f(s)dsj . (A7) 

The arbitrariness of uq can be exploited to find the solution of Eq. (A3) with the initial 
condition ip(k, 0) = ipo(k): 

V>(M)=exp|-« / u(k + [ f(s')ds'^\ ds\^ (k + [ f(s)ds). (A8) 



Appendix B: Statistics of the write-in variable 



The characteristic function Z(x) = J dx exp (ixx)V(x\pf) is 



Z(x) (x exp (iF X:T ) - 2A x X)T0 exp (iF x , T )lm.(A w ) - iy Xjtv exp (iF XjT )Re(A w ) 

v2 



(4,-o - Viu) ex P (iF x ,r) - o( d l F x,r ) exp (iF x , T 



B„ 



(Bl) 



with F Xjt = xx + [up{k + x/2) - Up{k - x/2)]i, z X|t = 9 x F X)t , y X)4 = d k F x ^ t v = (r - r )/2. 



The normalization is given by A = "P p0 st/cto,o- We note that exp (iF XtT ) is the characteristic 
function for x that the probe would have after a time r, had it not interacted with the 
system. 
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Appendix C: Measurement of a spin 1/2 



1. Exact expressions 



We use units of A&, kj A k — > k, A/A^ — > A, 6Ak — > 0. We assume the free Hamiltonian 
for the probe to be Hp = h 2 k 2 /2M. This position defines the dynamical scale: k 2 H = M/Htq, 
with r = f^clt dsg(s). 



Po(k — Act, k — Act) (<j\ pj \a) (a\ pi \a) 



+ e 2taXk/k xp (k — Act, k + Act) (-a\ p f \a) (a\ Pi \-a) 



(CI) 



For definiteness po(k, k') oc exp {—[(A; + fc') 2 /8 — (k — k r ) 2 /2nl]}, with n k < 2 the coherence 
scale. We let pi = (1/2)(1 + m • cr), pf = (1/2)(1 + n • cr), A = a • cr, with a a unit vector 
and |m| < 1, |n| < 1. The probability of a successful postselection is then 

1 



V, 



post 



1 + mana + e 2A2/K fc 2A2/fc » (m-n - m-an-a) 



(C2) 



The joint probability is 



1 



V{k,p f ) = - { (1 + m-a)(l + n-a)P (fc - A) + (1 - m-a)(l - n-a)P (fc + A) 



i '2 Afc 

. mn — mana) cos —5— 

k z 



H 



. ,'2Afc 
(mxnj-a sm( -yy- 



Po(fc) • (C3) 



The characteristic function is 



Z{6) = |^-<! (1 + m an a) cos (A0) + i (ma + n-a) sin (\6) 



'Ppost 2 



+ e ~ 2A2 K~ 2 +^ 4 ) 



m n — m an a cosh — tt- — 1 mx n -a smh 



k 2 



k 2 



(C4) 



with Z o (0) = exp [-fl 2 /2]. 



2. Coherent oscillations 



Equation (C3) shows that the off-diagonal terms of po(k,k') provide oscillating terms, 
with a period k osc = nk 2 H /\. If k osc exceeds the scale A& over which Vo(k) decays, the 
oscillations can not be discerned. If instead k osc <C A&, coherent oscillations are observed. 
In this case, we can not expand in A the exponential of the Hamiltonian phases exp (iT a ), 
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but should retain the full terms. This can be implemented through the transformation 
Pf — > UkPfUl = Pf(k), with the unitary operator 

U h = exp {i\u' P (k)T A - i^(fc)(r - r^i 2 }. (C5) 

For a spin 1/2, £4 is but a rotation around the direction of A. Accordingly, the weak value 
becomes a periodic function of k 

Tr{ Pf (k)A Pl } Tr{ Pf (k)A Pl A} , . 

Furthermore x t — >■ x and x t \k — > x in Eqs. (10)-(14). Once these prescriptions are applied, 
the probability in Eq. (13) approximates excellently the exact expression. 
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